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We theoretically analyze quasi-one-dimensional Bose-Einstein condensates under the influence of 
a harmonic trap and a narrow potential defect that moves through the atomic cloud. Performing 
simulations on the mean field level, we explore a robust mechanism in which a single dark soliton is 
nucleated and immediately pinned by the moving defect, making it possible to drag it to a desired 
position and release it there. We argue on a perturbative level that a defect potential which is 
attractive to the atoms is suitable for holding and moving dark solitons. The soliton generation 
protocol is investigated over a wide range of model parameters and its success is systematically 
quantified by a suitable fidelity measure, demonstrating its robustness against parameter variations, 
but also the need for tight focusing of the defect potential. Holding the soliton at a stationary 
defect for long times may give rise to dynamical instabilities, whose origin we explore within a 
Bogoliubov-de Gennes linearization analysis. We show that iterating the generation process with 
multiple defects offers a perspective for initializing multiple soliton dynamics with freely chosen 
initial conditions. 

PACS numbers: 03.75.Lm,67.85.De,67.85.Hj 


I. INTRODUCTION 

Solitonic wave excitations that maintain their shape 
during propagation are found in a large variety of phys¬ 
ical systems, ranging from hydrodynamics to modern 
telecommunication systems and even biological molecules 
[T]. The enormous technological advance of recent years 
has made it possible to prepare and observe solitons in 
ultracold atom experiments, primarily (but not only, see 
e.g. 0) in condensed bosonic ensembles near zero tem¬ 
perature mm- On the mean-field level, a Bose-Einstein 
condensate (BEG) is described by the Gross-Pitaevskii 
equation (GPE) [3], a nonlinear Schrodinger equation 
with cubic nonlinearity induced by the interatomic in¬ 
teraction. In one spatial dimension (ID), this equation 
is well-known to feature dark (bright) solitons for defo- 
cusing (focusing) nonlinearities, respectively [U |S] . Ex¬ 
perimentally, a highly elongated quasi-lD regime can be 
reached by tightly confining the atoms in the radial di¬ 
rection, effectively freezing out the transverse dynam¬ 
ics. Early experiments succeeded in preparing bright 
[SlIZl and dark mis] matter-wave solitons. Dark soli¬ 
tons (which are in the focus of this work) are charac¬ 
terized by a localized density minimum across which the 
phase changes by tt. These can be created by manipu¬ 
lating the condensate phase [SHm or density m using 
external potentials, see also the discussion in They 
can also form in the wake of a repulsive barrier dragged 
through the condensate [SHin] or in collisions of initially 
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separated atomic clouds m- Still, controllably creating 
a dark soliton at a desired position in a Bose-Einstein 
condensate is a challenging task in experiments. 

Once a soliton has formed, its dynamics can again be in¬ 
fluenced with external potentials. A dark soliton in a har¬ 
monically trapped BEG performs particle-like harmonic 
oscillations around the trap center [ISIIIIIIMIIIT]. The 
emission of sound waves due to acceleration of a soliton 
has been studied in [IHHII]. Thinking of the external 
force as a handle for controlling the soliton dynamics, 
a direct manipulation of the soliton motion with nar¬ 
row potential defects has been suggested in [^, where 
the interaction of a dark soliton with a pointlike impu¬ 
rity was analyzed within perturbation theory (see also 
and the corresponding studies for bright [55] and 
dark-bright solitons m)- Moreover, the possibility of 
dragging along dark solitons in a moving optical lattice 
potential has been demonstrated in |2S||5H|. In a simi¬ 
lar spirit, dragging of bright solitons in a discrete lattice 
model has been discussed recently |30j , while pinning and 
transporting quantum vortices with focused external po¬ 
tentials has been shown in [niisij. 

Here, we describe a method for controllably creating, 
dragging and releasing a dark soliton in a repulsively in¬ 
teracting Bose-Einstein condensate employing a tightly 
focused red-detuned laser beam (acting as an attractive 
potential for the atoms via the dipole force). Specifically, 
we study a trapped quasi-ID BEG under the influence of 
a moving Gaussian potential defect of attractive sign. 
When entering the atomic cloud from outside with this 
defect, a dark soliton can be created and at the same time 
pinned by the defect, such that it can subsequently be 
placed and released at an arbitrary position. Motivated 
by an analysis of the instantaneous energy levels in a sin- 
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gle particle model, a similar scheme for exciting nodes 
in the condensate wave function by traversing it with an 
attractive defect has been proposed in [231 [M]. A related 
protocol for the extraction of bright solitons from an at¬ 
tractive BEC has been suggested in |3S] (see also [25]b 
but there the defect was not assumed to act as an ex¬ 
ternal potential, but instead to cause a local variation of 
the effective atomic interaction through the mechanism of 
optical Feshbach resonance m- In our work, the defect 
acts as a single particle potential for the atoms. Dynam¬ 
ically manipulating (e.g. splitting) the entire BEC cloud 
with such optical “tweezers” is nowadays well established 
[3Sll2n]. If, instead, the light is focused to the compara¬ 
bly short length scale of a dark soliton, the same type of 
technology can be employed to manipulate the dynamics 
of a localized solitonic excitation. 

Our presentation is structured as follows. In Sec. |lT| we 
introduce the theoretical mean-field framework our study 
is based on and give a discussion of the results from a 
perturbative treatment of the potential defect, revealing 
in particular that a potential that is attractive to the 
individual atoms is effectively also attractive to a dark 
soliton. In Sec lIIIl we show results from numerical simu¬ 
lations, demonstrating the robust creation, dragging and 
release of a dark soliton. We quantify the fidelity of the 
creation process and vary the defect parameters to ex¬ 
plore the robustness of this protocol. If the dark soliton 
is pinned to the defect for long times, a dynamical in¬ 
stability may occur which we address in Sec |IV| making 
the connection to a Bogoliubov-de Gennes linearization 
around the corresponding stationary solution of the GPE. 
We briefly conclude and point to further perspectives in 
SecjV]. Details of the perturbation theory are defered to 
appendix [^ 


II. SETUP AND RESULTS FROM 
PERTURBATION THEORY 

We investigate a quasi-ID Bose-Einstein condensate of 
a single atomic species with repulsive interaction at zero 
temperature, described by the Gross-Pitaevskii equation 
[3] . We assume strong harmonic confinement in two spa¬ 
tial directions, taking the trapping potential to be of the 

form Vi 3 D(i', t) = ^’^|^(j/^-l-z^)-l-P(a;,t), where m denotes 
the atomic mass, uj± the frequency of the transverse os¬ 
cillator potential and V{x,t) models a potential in the 
longitudinal direction. Assuming that the transverse dy¬ 
namics is fully frozen out and thus transversally the con¬ 
densate wave function remains in the oscillator ground 
state, one can integrate out the y- and z-directions and 
is left with the effectively ID Gross-Pitaevskii equation 
01 

+ V{x,t)tp + gmlV'lV = (1) 

where tp = tp{x,t) denotes the longitudinal part of the 
wave function and the nonlinearity coefficient giu = 


2ahw± with the s-wave scattering length a > 0 [3j . Mea¬ 
suring length, time, energy and density {ip^x, t)\'^ in units 
of ax = i/S/mwx, fiuj± and (2a)“^, respectively, 
Eq. Q is cast into the dimensionless form 

-^dl'ip + V{x,t)tj;+\tp\'^ilj = idt'ip ( 2 ) 


which we will work with in the following. 

The corresponding stationary equation is obtained by 
factorizing ipi^x, t) = exp(—i/it) with g the chemi¬ 

cal potential. Our focus here will be on a longitudinal 
potential that consists of a static harmonic part (whose 
frequency wy is small compared to that of the transverse 
trap) plus a Gaussian of fixed height and width, but mov¬ 
ing in time, i.e. in dimensionless units 


V{x,t) 
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Here, Xa{t) specifies the trajectory of the Gaussian im¬ 
purity, while a and Vq set its width and amplitude, re¬ 
spectively. For convenience, the aspect ratio will be fixed 
to = Wy/w^ = 0.04 in the following. Our results can 
be transfered to other aspect ratios by a straightforward 
rescaling of Eqs. (HI- We simulate the time evolution 
of the condensate by propagating Eq. © with a fourth- 
order Runge-Kutta integrator. The initial state is chosen 
to be the ground state at a given y, obtained from the 
stationary ID-GPE by an adapted Newton method m- 
The main objective of this work is the controlled forma¬ 
tion and dragging of a dark soliton. By analogy with 
previous work on vortices in 2D |32j . and the somewhat 
intuitive picture that a density dip such as the dark soli¬ 
ton may offer the possibility of pinning it with a narrow 
repulsive barrier in its center, one may conjecture that 
Vo > 0 will be the favorable parameter regime for our 
purposes. It turns out, however, that this is not the 
case and an attractive Gaussian impurity with Vq < 0 is 
much more adequate to drag along the soliton. This can 
be seen on the level of dark soliton perturbation theory 
m. resulting in an approximate particle-like equation of 
motion for the soliton center in the presence of a weak 
external potential. This has been worked out for a dark 
soliton perturbed by a Ashaped impurity potential in 
m, with the result that the impurity is attractive (re¬ 
pulsive) to the soliton if it is attractive (repulsive) to the 
atoms in the condensate. We have extended this analysis 
to our impurities of Gaussian shape (see appendix [A| for 
the details) and And that the overall result persists: Dark 
solitons are effectively attracted by a Gaussian of Vq < 0. 
Specifically, the soliton center xq follows the equation of 
motion 




dW 

dcro ’ 


( 4 ) 


where the effective potential W (xq) is predicted from per¬ 
turbation theory, see appendix [A] for the details. Fig. 
shows the resulting effective potentials for varying widths 
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X (in units of Oj^) 

Figure 1. (Color online) Effective potentials W(xo) for the 
soliton position coordinate xq as obtained from soliton per¬ 
turbation theory. In each case, the underlying atomic poten¬ 
tial consists of a harmonic trap and a Gaussian of variable 
width a and fixed attractive amplitude Vb = — 1 centered at 
a: = 0 (black circles for repulsive potential Vb = 1), /a = 1 
throughout. 
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of the Gaussian impurity (located at a: = 0), keeping the 
amplitude fixed. Clearly, a potential minimum for the 
soliton dynamics is found for all cases with Vb < 0, while 
it turns into an unstable maximum for Vq > 0 (black cir¬ 
cles in Fig. [^. In the vicinity of the fixed point x = 0, 
the effective potential is strongly shaped by the Gaus¬ 
sian, while further away from the center it asymptotes to 
n^x^/4, half the bare potential of the trap, yielding the 
characteristic soliton oscillation at 17/>72 as expected in 
the absence of the Gaussian 0117 ]. 

This perturbative treatment suggests that for dragging 
along a dark soliton with a Gaussian impurity, one should 
choose the impurity as attractive for the atoms. Conse¬ 
quently, we will focus on Vo < 0 in the following. It 
should be noted that the parameters used in the rest of 
this work are mostly out of the range of validity of the 
soliton perturbation theory. 


III. GENERATING, PINNING AND 
DRAGGING A DARK SOLITON 

In this section, we demonstrate the possibility to gener¬ 
ate, pin and drag along a dark soliton with an attractive 
Gaussian impurity entering the BEG cloud from outside. 
An example of this is shown in Fig. da), displaying the 
spatio-temporal evolution of the atomic density under the 
influence of the impurity potential. The dimensionless 
chemical potential is chosen as /i = 1, corresponding for 
instance to a condensate of around 3300 sodium atoms 
under a transverse confinement of u}± = 2 tt x 200 Hz. 
For these parameters, the resulting healing length in the 
center of the cloud is close to one micron. The white line 
indicates the trajectory of the Gaussian that moves lin¬ 
early into the BEG cloud towards x = 1. After staying 


Figure 2. (Color online) Generating, dragging and releasing 
a dark soliton. The white line indicates the trajectory of 
the Gaussian impurity (parameters Vo = —12 ,(t = 0.1). At 
t = 129 the impurity is switched off. (a) Density \ip{x,t)\'^. 
(b) Snapshots of the phase angle profile at different times. 


stationary at this point for a time interval of At = 10, it 
is switched off. The generated soliton can be identified 
already at an early stage. When the impurity enters the 
cloud, the characteristic density minimum as well as the 
phase shift close to tt (see Fig. IJb)) are created almost 
immediately. The soliton follows the motion of the impu¬ 
rity and is dragged along towards x = 1 where it is held 
for At = 10. When the Gaussian potential is switched 
off, the soliton is released and starts to oscillate in the 
harmonic trap. On a perfect Thomas-Fermi background, 
the frequency of this solitonic oscillation is expected to 
be r7/-\/2 0[T7]. In the present simulation, the Gaussian 
impurity also slightly excites the collective dipole mode, 
causing a center-of-mass oscillation of the entire cloud at 
the trap frequency O. A two-sine fit to the soliton tra¬ 
jectory reveals a superposition between the particle-like 
soliton oscillation at r7/-\/2 and the collective dipole os¬ 
cillation at n, validating the dark soliton character of the 
created excitation. 


We have confirmed in our simulations that the gener¬ 
ation, immediate pinning and dragging of a dark soliton 
as shown in Fig. is successful in a wide range of param¬ 
eters. Specifically, while the velocity of the impurity is 
relevant (see also below), the exact trajectory is not, and 
the process works equally well with curved trajectories 
of the impurity and shorter or longer hold times prior to 


release (see, however, the discussion in Sec. IVI. 
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eration process, we need to identify a scenario in which 
a precise definition of the ideal desired outcome can be 
given. While on a homogeneous background there is a 
notion of a perfect grey (moving) soliton, an extension of 
this accounting for the inhomogeneous density induced 
by the trap is available on an approximate level only. In 
contrast, the profile of a fully stationary (black) soliton 
in the trap can be obtained unambiguously by solving 
the stationary GPE. Moreover, it can be expected that 
at parameters, at which an undisturbed initialization of 
a black soliton is possible, an off-center release from the 
impurity and the subsequent acceleration due to the trap 
(similar as seen in Fig. will produce a clean grey soli¬ 
ton. Thus, we proceed to quantitatively evaluate the 
success of the soliton generation process by focusing on 
the preparation of a black soliton in the trap center. To 
do so, we choose Gaussians of different parameters that 
enter the cloud at a given velocity and move towards 
X = 0, are held there for a while and then turned off. For 
each of these runs, we compare the resulting final state 
to that of a stationary black soliton, our target state. 
This target state (/>bs( 3^) (with the same squared norm 
N as the wave function in the simulation) is computed 
separately by numerically solving the time-independent 
GPE. Then, for any time t after the potential has been 
turned off, we can calculate the overlap 


S{t) = 


1 

iV2 


2 


dx (l)Qs{x)tp{x,t) 


(5) 


This quantity is then averaged over a time-interval 
[trjtf], ranging from the release time tr (when the im¬ 
purity is switched off) to a final time tf\ 

S= ^ f dtS{t). (6) 

Jtr 


By the Gauchy-Schwarz inequality, 0 < S' < 1. If the 
dynamics caused by the impurity results in a perfect sta¬ 
tionary black soliton, one has |V'(x,t)| = |(/>Bs(a;)| for all 
t > tr (taking advantage of the fact that the target state 
is stationary), and thus S = 1 would correspond to a 
perfect fidelity of the creation process. Smaller devia¬ 
tions from this indicate dynamics of the soliton and/or 
the bulk of the cloud after the impurity has been turned 
off, while S ^ 1 suggests that the generation of a single 
dark soliton has completely failed. Fig. shows results 
for the fidelity S' as a function of width tr and height 
Vb < 0 of the Gaussian impurity. The 1/e^ width of the 
Gaussian is given hy w = 4a and ranges from 0.24 to 
1.2 here. This is to be compared to the healing length in 
the center of the cloud given by ^ « 0.7. The impurity 
moves towards the trap center on a linear trajectory, as 
in Fig. 1^ at a velocity that is fixed to u = 0.0925 in this 
set of simulations (measured in units of a±/uj±}. It is 
then suddenly stopped and held at x = 0 for At = 10, 
before being switched off. The quantity S is then ob¬ 
tained by averaging S(t) from the subsequent dynamics 
over an interval tf — tr = 481. 



-20 -15 -10 


Vq (units of 

Figure 3. Black soliton generation fidelity 5" as a function 
of the potential strength Vb and width a. The Gaussian en¬ 
ters the cloud at a velocity v = 0.0925, the initial chemical 
potential /r = 1. 


Most notably, there is an extended parameter region of 
substantial fidelity 5 > 0.9. A small or intermediate po¬ 
tential width a together with a large or intermediate po¬ 
tential strength |Vb| is applicable for the controlled gen¬ 
eration, dragging and holding of a single dark soliton. 
An exemplary plot at the same parameter values as in 
Fig. i (but now with the final position of the Gaussian 
at X = 0) is shown in Fig. ^a) for the parameter set 
(Vo = —12,cr = 0.1). For this comparably narrow im¬ 
purity, the final state is close to a stationary black soli¬ 
ton, but both the dipole mode and the soliton oscillation 
mode are slightly excited. In contrast, turning to the 
parameters (Vb = —7,cr = 0.28), the fidelity is roughly 
the same, S ~ 0.94, but the deviations from the black 
soliton state are of a different kind, see Fig. Qb). Here, 
the soliton itself is closer to stationary than in Fig. ia), 
but the background is excited at higher frequency modes 
and more disturbed by density waves. The comparison of 
Figs. I^a) and (b) illustrates that the fidelity S is sensi¬ 
tive to different types of remaining excitations around the 
target state (both particle-like oscillations of the soliton 
and collective oscillations in the bulk), and that one has 
some freedom in reducing either the particle-type soli¬ 
ton dynamics or the background excitations by tuning 
the parameters of the Gaussian. In both regions of pa¬ 
rameter space, the dark soliton characteristics of the in¬ 
duced density minimum are clearly observed; specifically, 
we have checked the oscillation at n/-\/2 for off-centered 
release. Finally, Fig.Qc) shows a simulation at parame¬ 
ters (Vb = —8, cr = 0.19) that yields a particularly large 
fidelity of S' « 0.98. Here, the evolution is similar to that 
shown in Fig. Sb). but the background excitations are 
further suppressed. 
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Figure 4. (Color online) Density evolution \'ip{x,t)\'^ for some 
of the simulations underlying Fig. The parameters of the 
Gaussians are {Vo, a) = (—12,0.1) (a), (—7,0.28) (b) and 
(—8,0.19) (c), respectively. 


So far, we have not addressed the role of the ve¬ 
locity at which the Gaussian defect is moved through 
the condensate. For the two parameter sets {Vo,a) = 
(—12,0.1), (—7,0.28), we have performed simulations for 
a range of defect velocities v. The resulting fidelity S 
for varying v is shown in Fig. For both parame¬ 
ter sets, we find a robust plateau of large S at small 
V. In particular for the deep and narrow Gaussian with 
(Vo, cr) = (—12,0.1) the fidelity is close to 1 for u < 0.08. 
Increasing the velocity, for this parameter set a sharp 
drop of S is observed at u « 0.16. Beyond this velocity, in 
addition to a rather strong excitation of the background, 
further dark solitons are generated of which at most one 
is pinned by the Gaussian, which drastically reduces the 
overlap with the single black soliton state. This param¬ 
eter regime could be of interest in its own right when 
moving the focus towards multiple soliton physics, but 
the high degree of control over the single soliton creation 


Figure 5. (Color online) Black soliton creation fidelity as a 
function of the impurity velocity for two different sets of pa¬ 
rameters of the Gaussian potential. 


process is lost there. In contrast, a comparable critical 
drop in the fidelity is not observed for the Gaussian with 
{Vo, a) = (—7,0.28), where S remains relatively large for 
an extended range of v. For these parameters, going to 
even larger velocities than shown in Fig. we observe a 
trend of decreasing S, but caused by enhanced excitation 
of the background and the soliton oscillation mode, in¬ 
stead of multiple soliton formation. 

In this work, we do not aim to explore the full velocity de¬ 
pendence and the dynamical details of the nucleation pro¬ 
cess in which the dark soliton is created when the defect 
enters the cloud. Generally, the formation of excitations 
in a superfluid simultaneously exposed to a moving de¬ 
fect and a trapping potential (which leads to an inhomo¬ 
geneous background density) is an interesting and timely 
subject to study in its own right, see for instance the 
recent experiment on vortex shedding in quasi-2D con¬ 
densates m- Even on a homogeneous background, pre¬ 
dicting the critical velocity above which a defect causes 
the creation of nonlinear excitations is an intricate prob¬ 
lem (see for instance the extensive discussion of vortex 
nucleation in m) that is subject to ongoing research 
gam]. Some related results on soliton formation in ID 
(mostly focusing on repulsive defects) are available, e.g. 
[nmaisiiis], but in our setup additional complications 
due to the inhomogeneous background density traversed 
by the defect, the attractive sign and the comparably soft 
boundary of the Gaussian potential gT] and possibly also 
the ID reduction m will require separate investigations 
that are beyond the scope of the present study. 


IV. INSTABILITIES FOR LONG HOLD TIMES 

In the above simulations, the soliton was created and 
dragged by the moving impurity, then placed at a de¬ 
sired position in the trap and held there for a compa¬ 
rably short time {At = 10) before being released. Sub¬ 
stantially extending this hold time reveals an additional 
effect. For certain parameter values we observe instabil¬ 
ity phenomena in the dynamics of the pinned soliton, 
see Fig. We find that the soliton performs micro- 
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Figure 6. (Color online) Controlled creation of a dark soliton, 
followed by a long hold time. The parameters of the Gaussian 
are Vb = —12, a = 0.08, v = 0.0925 (before reaching x = 0). 
(a) Density evolution \'tp{x,t)\^. The piece-wise linear trajec¬ 
tory of the impurity is marked by a white line. The white 
box highlights the first region of dynamical instability, (b) 
Position of the density minimum (as a measure of the soliton 
center) as a function of time in the time interval marked by 
the box in (a). 


oscillations around the Gaussian potential during the 
hold interval (similar oscillations are observed in all our 
simulations). Here, however, at t « 140 the amplitude 
of the micro-oscillations starts to increase strongly, be¬ 
fore it decays again after t 210. A similar increase 
and decrease is observed again at around t > 300. We 
conclude from this that long hold times may give rise 
to undesired effects when aiming for a stationary soli¬ 
ton. Inadvertently releasing the soliton during a period of 
enhanced micro-oscillation will yield a comparably large 
momentum of its particle-like motion and correspond¬ 
ingly a large amplitude of its subsequent oscillations in 
the trap. Thus, the hold interval may crucially affect the 
final state of the soliton preparation process. Compa¬ 
rable dynamical instabilities of a dark soliton under the 
influence of a narrow external potential have been related 
to sound emission caused by the repeated asymmetric de¬ 
formation of the oscillating soliton due to the impurity 
m. see also [milHlIlHM]. 

To obtain further insights into the dynamical insta¬ 
bilities due to the Gaussian potential, we employ a lin¬ 
earization analysis. Let us assume that the state that 
originates from the Gaussian entering the cloud and 
moving towards a: = 0 is close to the corresponding 
stationary black soliton state (now in the presence of 


the Gaussian). Then, information about its stability 
is encoded in the Bogoliubov-de Gennes (BdG) exci¬ 
tation spectrum, obtained by adding a small deviation 
5'ip{x,t) = to the station¬ 

ary dark soliton state, linearizing the GPE in Sip and 
solving the ensuing eigenvalue problem for u [1]. Here, 
/r denotes the chemical potential of the stationary solu¬ 
tion. Frequencies having nonzero imaginary part indicate 
an instability of the intial state as they induce exponen¬ 
tial growth of a generic small perturbation. Such com¬ 
plex modes may emerge from collisions of normal and 
anomalous modes (modes with positive or negative en- 
ergy/Krein signature, respectively) [H |52]. Here, the 
BdG spectrum of the dark soliton state is expected to 
exhibit a single anomalous mode that is related to its 
particle-like motion [^. If this becomes resonant with 
one of the background modes as a parameter is tuned 
(for instance the width of the Gaussian), this may lead 
to instability of the state. Indeed, we observe these ef¬ 
fects in the BdG spectrum of the black soliton state with 
a Gaussian potential placed in its center. We fix the 
norm of the wave function and the amplitude Vq = —12 
of the Gaussian to the same values as in Fig. |^and scan 
a. The resulting spectrum as a function of a is shown 
in Fig. I^a) (by the Hamiltonian symmetry, if w is in the 
BdG spectrum, then so are — w, w* and —w*, so we can re¬ 
strict to positive real and imaginary parts in the figure). 
Here, even for moderately small a the large amplitude 
of the Gaussian strongly shifts the anomalous solitonic 
oscillation mode away from its value H/-\/2 expected in 
the harmonic trap only. Increasing a leads to a further 
increase of the anomalous mode frequency, causing sub¬ 
sequent collisions with background modes that result in 
complex quartets, signalling oscillatory instability. The 
width a — 0.08, as used in the simulation of Fig. in¬ 
deed lies at the edge of such a region of instability. Let us 
at this point return to the perturbative regime of small 
|Vol that was discussed in Sec. |h] Numerically calculat¬ 
ing the BdG spectrum as a function of a also reveals the 
emergence of complex instability bubbles in this regime, 
see Fig. m for an example. The solitonic perturbation 
theory does not account for background excitation modes 
and is not capable of predicting these instabilities. Inter¬ 
estingly, however, the linearization frequencies predicted 
from the perturbation theory (by linearizing the effective 
potential W{xq) around xq = 0) quite accurately cap¬ 
ture the real part of the unstable BdG modes even in the 
instability regions. In regions of stability, the frequency 
from the perturbative approach is close to the anoma¬ 
lous mode in the BdG spectrum, as expected. As is to 
be expected, the agreement slightly deteriorates for large 
a, where the overall perturbation due to the Gaussian 
effectively becomes stronger. 
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Figure 7. (Color online) BdG linearization spectra of the 
dark soliton state in the presence of a Gaussian potential in 
its center. All BdG frequencies are given in units of Blue 
crosses, red circles and green asterisks denote the normal, 
anomalous and complex modes, respectively. The chemical 
potential is fixed at ^ = 1.15, yielding the same norm as the 
ground state at ^ = 1. (a) Vb = — 12, revealing a complex 
mode at (T = 0.08, corresponding to the parameter set of 
Fig. (b) Vo = —0.25, including also the linearization fre¬ 
quency predicted from the soliton perturbation theory (black 
dots). 


V. PERSPECTIVES AND CONCLUSIONS 


We theoretically investigated the possibility to control- 
lably generate, drag, hold and release a dark soliton in 
a quasi-ID Bose-Einstein condensate using a Gaussian- 
shaped impurity as could be implemented with a focused 
laser beam. The time-dependent Gross-Pitaevskii equa¬ 
tion containing the trap potential as well as the Gaussian 
impurity was propagated in time to obtain the spatio- 
temporal evolution of the condensate wave function when 
disturbed by the moving defect. On a perturbative level, 
we found that if the Gaussian is attractive (repulsive) to 
the atoms, then it is effectively attractive (repulsive) to 
the dark soliton as well, thus suggesting the use of an at¬ 
tractive impurity (red-detuned focused laser) for holding 
and dragging dark solitons, in contrast to the pinning of 
vortices at repulsive barriers. 

We demonstrated that by entering the atomic cloud with 
an attractive Gaussian one can create a single dark soli¬ 
ton that immediately sticks to the defect and can be 
controllably placed and released at a desired position in 
the condensate, showing the expected characteristics of 
a dark soliton after release. Detailed investigations re¬ 
vealed an extended range of model parameters (such as 
the width and amplitude of the Gaussian) for which this 
mechanism is successful, thereby underlining its robust 



100 200 300 400 500 600 

t (units of wj^) 


Figure 8. (Golor online) Generating two solitons with two 
impurity potentials moving independently on the trajectories 
indicated by white lines. Golors encode the density \'tp{x,t)\^. 
For both Gaussian potentials, Vb = —14, a = 0.09. 


nature. As a drawback, the width of the Gaussian must 
be relatively small, comparable in size to the soliton heal¬ 
ing length, requiring a much tighter focus than in previ¬ 
ous experiments m- For instance, in the case of ^^Na 
atoms and a transverse confinement of a;_L = 27rx 200 Hz, 
a dimensionless value of ct = 0.2 in our simulations (for 
which we observe particularly successful soliton genera¬ 
tion and control) translates into a 1/e^ beam width of 
w = 4cro_L « 1.2/rm, close to the central healing length, 
while the wave function norm in our simulations trans¬ 
lates into a relatively small number of around 3300 atoms. 
This is a challenging requirement, but not out of reach, 
given that optical systems with sub-micron resolution are 
already employed in present-day cold atom experiments 

[53 EH- 

Moreover, our studies suggest that long stationary hold 
times of the defect are not favorable for the controlled 
generation of black solitons, due to the possibility of dy¬ 
namical instabilities that may lead to a spontaneous in¬ 
crease of the micro-oscillation amplitude of the soliton 
around the Gaussian. This was related to correspond¬ 
ing complex modes arising in the Bogoliubov-de Gennes 
spectra, and linked to the linearization results from the 
perturbative approach. 

The protocol for simultaneous generation and holding of 
the dark soliton as described here is appealing, since cap¬ 
turing an existing soliton in a BEG cloud would be a 
much more difficult task. Even if a suitable pinning po¬ 
tential is available, catching the soliton requires informa¬ 
tion about its time-dependent position, which is hard to 
obtain given the destructive measurement schemes. Fur¬ 
thermore, we point out that if more than one laser beam 
is available, the soliton creation scheme described herein 
can immediately be cascaded to generate two or more 
dark solitons at predefined positions as shown in Fig. [^ 
cf. also [33]. The two dark solitons are created and 
trapped by their respective defect potential, and when 
released they exhibit the expected particle-like collision 
dynamics, cf. USES]. 

In this work, we have not addressed the details of 












































the soliton nucleation process occuring in the low den¬ 
sity wings of the cloud and its dependence on the defect 
velocity. This aspect of the problem promises to be an 
interesting topic for future studies. Also, it would be 
desirable to perform simulations of the soliton creation 
and dragging protocol in the framework of the full three- 
dimensional Gross-Pitaevskii equation, checking for pos¬ 
sible transverse excitation effects that are not captured 
within the dimensional reduction. 
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Appendix A: Soliton perturbation theory 


In this appendix, we outline the derivation of the effec¬ 
tive potential for the soliton in the harmonic trap, per¬ 
turbed by a Gaussian impurity. We follow the presenta¬ 
tion in EU, but generalize the Dirac (^-potential used in 
that work to a Gaussian as in Eq. ^ . By the same argu¬ 
ments as in E2], the Thomas-Fermi-like background pro¬ 
file Ub{x) of the condensate ground state is approximated 
by Ub{x) = uo+ftra.p{x)+fg{x), where uq is the maximum 
background amplitude {uq = in the Thomas-Fermi 
limit considered here), /trap(a^) accounts for the modified 
shape due to the harmonic trap and fg{x) incorporates 
the perturbation by the Gaussian impurity. Explicitly, 


/trap(^) — p, htrap(^) with Ftrap — ^ X /2, as in 
ZUq 


and 



(Al) 

Here, we can recover the result of E2] by taking the Dirac 
limit Vo = bl{\/^a) and cr —?> 0 for a fixed b, resulting 
in fg{x) = -{b/2)exp{-2uo\x\). 

With the generalized Ub{x), we follow the further steps 
in El]. The dynamics of the dark soliton on top of the 
Thomas-Fermi-like background is investigated with the 
Ansatz 


'ijj{x,t) = Ub{x)e 


(A2) 


where u(x, t) represents a dark soliton on this back¬ 
ground. Inserting Eq. (A2) into the time-dependent GPE 


leads to a perturbed defocusing nonlinear Schrodinger 
equation for the soliton function v, with a perturbation 
term that depends on /trap and fg. Making an Ansatz for 
V in the form of a dark soliton solution of the defocusing 
nonlinear Schrodinger equation, but with its position and 
phase angle slowly varying in time, one can employ the 
adiabatic perturbation theory for dark solitons of E3| to 
obtain the desired equation of motion of the soliton cen¬ 
ter Xo(t), which in our case reads as 


d^xp 

dt^ 


1 dVtrap 

2 dx 


-k 




dW 
dxp ’ 


(A3) 


where the integrands 


Ei(x) 

F2{x) 


= I — 1 — erf 


O' / 


— l^-2uo+^ 


^-2uox 


X {tanh[uo(x — xp)] — 1} sech^ [mp(x — xp)] 


= — 1 -k erf 


V2 


(2up + 


c^UoX 


X {tanh[up(x — Xp)] -k 1} sech^ ['ap(x — Xp)] 


From Eq. |A3| we can numerically compute the effective 
potential VF(xp) (as shown in Fig.[^ and the linearization 
frequency around its fixed point at xp = 0. 
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